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1. INTRODUCTION 

Transformers are a vital and expensive component of electrical power systems. The protection of 
transformers is essential to achieve high quality performance in modern electric networks. The winding faults 
on power transformers caused by high sensitivity, selectivity and fast response will be deterred by the 
differential protection relay [1]. Transformers with a high capacity are usually protected by a harmonic 
restraint percentage differential relay [2]. The faults that occur in transformers are either: internal incipient 
faults or internal short circuit faults [3]. Short circuit faults constitute about 70-80% of transformer failures, 
and these can be (1) phase-to-ground fault or (2) turn-to-turn faults [4]. 

Due to the nonlinearities in the transformer core, or in the CT core or in both case a substantial 
differential current may flow, when there is no fault. These false differential currents are generally sufficient 
to cause a percentage differential relay to trip. Some of these phenomena are magnetizing inrush current 
during energization or fault removal, transformer over-excitation and CT saturation. However, in modern 
digital relays some algorithms were developed to avoid the false tripping in percentage differential relays. 
Some of these algorithms have been developed based on different techniques, such as, artificial neural 
networks, fuzzy logic, wavelet transforms and principal component analysis [5-8]. These approaches have 
limitations that may affect their speed of operation, dependability and security. They also require complex 
computations and are susceptible to changes in transformer parameters. 


Journal homepage: http://ijpeds.iaescore.com 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 O 1089 


The differential protection on magnetizing inrush currents can be averted by many effective 
methods [9-12]. One of those methods is making the phase voltage as a control signal. And this way 
guarantee that the protected transformer is explicitly decreased in the terminal voltage. Also, implementing a 
core-flux or transformer which is equivalent circuit for electing inrush from internal fault currents could be 
another remarkable method. Those ways could achieve major success concerning their applications if the 
transformer terminal voltage measured and the measurement process will be quite costly [13]. Additionally, 
computational burden on the differential relay would be increased so that causes a slow motion in winding 
short circuits [14-16]. 

Phase angle difference technique, proposed by [1] is used to avoid limitations in previous 
techniques, and to improve the reliability and fast response in percentage differential relay. In this paper to 
discriminate magnetizing inrush currents and short circuit faults in power transformer, the PAD technique 
with percentage differential relays will be modeled and simulated using Simulink package in MATLAB 
software. Also, the real circuit model of power transformer and CTs transformer will be considered in the 
simulation model in both cases. 


2. METHODOLOGY OF DIFFERENTIAL RELAY WITH PAD 
2.1. Percentage differential (87T) characteristic 

Power transformers that are rated above 10 MVA are most commonly protected with percentage 
differential relay (87T) for the purposes of avoiding internal short circuits. The 87T-relay has been found to 
offer important ground and speed phase protection for 2 and 3 winding power transformers [17]. The relay 
features an additional harmonics restraint unit in which massive transformer magnetizing inrush is present as 
well as ratio matching taps. The unit prevents relay operation on transformer magnetizing inrush current. 
The second harmonics restraint unit is factory calibrated to restrain 15% second harmonic current, but may 
be adjusted if required. An unrestrained instantaneous unit which operates on magnitude of difference current 
is provided to back up the percentage differential unit. The unit is adjustable from approximately 8 to 20 
times tap [18]. The characteristic of percentage differential relay is providing following in Figure 1 [19, 20]. 
The characteristic In this case the differential current is false, they still can cause a trip in a percentage 
differential relay. Sometimes during energization inrush current is magnetized. The situation can also occur 
during transformer over-excitation, or fault removal, or CT saturation. However, recently, algorithms have 
been developed to avoid the false tripping in percentage differential relays. Various techniques such as 
principal component analysis, fuzzy logic and artificial neural networks helped to to develop these 
algorithms. These techniques have their parameters in Figure 1 are as follows: 
laip : is Magnitude of 50 Hz component for differential current. 
Tres: is magnitude of 50 Hz component for restraint current. 
Topi: is the minimum operation current. 
Iop2: is the adjusted minimum operating current 
Lees, min. is Minimum restraint current. 
k: gradient of the functional characteristic which are, 10, 20 or 40%. 
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Figure 1. Percentage differential relay (87T) characteristic 


Lag with transformation ratio of n is defined by (1), while Ies at the same ratio is given by (2). 
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laiff,o = lio = nh (1) 
Les, 7 g 5 nly g|/2 (2) 


Where, 
I, ø: is transformer primary current at phase Ø (namely, A, B or C). 
leg: is transformer secondary current at phase Ø (namely, A, B or ©). 


The percentage differential operating characteristic prevents operation until the differential current 
is greater than a certain percentage of through current. (3) as shown provides the operating characteristic of 
one phase differential element. The threshold value Iopl biases the differential operating current. 
The threshold value must be chosen based on the magnetizing current’s magnitude, and the differential 
current, which results from on-load tap-changing which occurs during normal transformer loading conditions 
[21, 22]. The second inequality models the slope of the transformer differential relay. 


laiff 2 lop 3) 
lai f = 4 Gee — Desmi) + lopi 


The minimum operating current Iop1 in amperes can be calculated with the following relation [11]: 


lpi =—XRXT (4) 


1+ 


R: is adjustable min. restraint setting, Re[1-3]. 
T: is tap setting in amperes. 


2.2. Fault discrimination techniques in power transformers 
2.2.1.Transformer current transient components technique 

The described technique is utilized to discern between external and internal faults within power 
transformers. The phase transient currents (all three) of the transformer go through conversion into the modal 
current components by way of Clarke's transformation to give ground mode I0, areal mode I1 as well as areal 
mode I2. High and low tension sides of the transformer give transient modal currents which are used with the 
Fault Discrimination Equation (FDE), which in turn is able to discern internal from external faults based on 
polarity of its output. The polarity of the FDE, for internal faults, will register negative whereas it will be 
largely positive for external faults. The following equation provides the fault discrimination equation: 


FDE = average{Uo)in* Uot + Uden © Ctr + Uden * Cader (5) 


Where, 

(Io)n , (Io)u : are the currents (transient ground mode) on high and on low tension sides. 
(IDn, U1)u : are the currents (transient aerial mode) 1 on high and on low tension sides. 
(D) , (l)u : are the currents (transient aerial mode) 2 on high and on low tension sides. 


2.2.2.Fault detection process in power transformers applying the Discrete Wavelet Transform 

The Discrete Wavelet Transform technique is used for on-line fault detection within power 
transformers. Wavelet technique is a time-scale domain approach, which is applied to locate short circuit 
faults and incipient faults through comparison of performance during normal operation of the power 
transformer. The following steps outline the fault detection algorithm [23]: 

Step 1: Obtain and record current signals from power transformer terminals. 

Step 2: Analyze the signal applying a wavelet taken from a wavelet family, for the required level of 

decomposition. 

Step 3: Assess the approximate coefficients and detail of Discrete Wavelet Transform with time 

through firstly plotting a sample—coefficient graph. 

Step 4: Locate the fault via wavelet coefficient interpretation. 

Step 5: Differentiate fault current (incipient), from internal short circuit current and ordinary 

operation current. 
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Discrete wavelet transform is also a technique that is used. It can be applied together with back 
propagation neural networks and is applied to classify different internal fault types in a three-phase 
transformer [24]. 

2.3. Phase angle difference (PAD) principle technique 

To differentiate the internal fault current from a different disturbance in a power transformer, 
the PAD in the fundamental frequency components will be used as the guide [25, 26]. It is not necessary that 
the exact PAD is used. The location of O relative to the PAD scheme operating characteristic can be 
estimated to find the status of a given power transformer. Figures 2 and 3 highlights that the operating mode 
of the power transformer is the main determinant of the principle of Phase Angle Difference. Under normal 
operation as well as in external fault, the line currents flow are parallel with the variance in phases being 
almost zero as shown Figures 2(b) and 2(e) respectively. However, during the magnetizing inrush currents, 
the PAD is approximately 90o since it is largely inductive as seen in Figure 2(a). Subjecting the transformer 
to winding short current reverses either, if the transformer is exposed to a winding short circuit, either I1 or 
nI2 current. In this case the PAD becomes greater that 90o in turn-to-turn fault and turn-to-ground fault as is 
visible in Figures 2(c) and 2(d). Now, through using Idiff and Ires, we can estimate PAD without placing 
much strain on the relay with great computations which would in turn slow its response. Figure 3 shows how 
to estimate PAD in a power transformer on different operation modes. 





Phase/deg. Phase/deg. Phase/deg. Phase/deg. Phase/deg. 


0.05 0.15 0.20 
(e) Time (s) 


Figure 2. PAD between the 50Hz components of phase A respective currents. (a) Magnetizing inrush. (b) 
Normal loading. (c) Turn-to- ground fault (d) Turn-to-turn fault and (e) external fault (AC-G) 
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Figure 3. Phasor diagram for computing PAD respective line currents: Normal loading, external faults, 
magnetizing inrush and internal faults 
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3. MODELING OF POWER SYSTEM COMPONENTS 
3.1. Simulink model of power system 

Figure 4 illustrates the power system’s single-line diagram used in the study. The sources FE, 0° 
and E,L @behind series impedances represent transmission and distribution networks. BA and BB are circuit 
breakers with resistance 0.001 Q. 


T:25MVA 
Bus1:138kV 138kV Y/13.8kV Y Bus2:13.8kV 


E.Z0° : L, / E LO’ 






Internal fault 


87T & PAD 


Figure 4. Single-line diagram of simulated power system 


External fault 





Figure 5 shows the Simulink model of the elements of the power system with the protection relay. 
The power system parameters in Figure 5 (a) are given in the Appendix. The power transformer in Simulink 
model in Figure 5(a), can be noted as in Y/Y together showing a grounded neutral together with a 25MVA 
rated capacity and 138kV/13.8kV rated voltages. As shown in Figure 5(b), the power transformer is modeled 
for calculating different operating modes, such as: internal faults, normal operation and magnetizing inrush 
currents. We can simulate the operations through power transformer control switches S1, S2, S3, SOI and 
S02 with BA and BB circuit breakers. 
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Figure 5. Power system with 87T relay and PAD scheme: (a) three phase diagram (Simulink model) and (b) 
Power transformer model (phase A) 
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The cores and windings of transformer is modeled by equivalent circuit illustrated in Figure 6(a). 
The (current— flux) characteristics of power transformer and CTs are shown in Figures 6(b) and 6(c). 
The calculations of transformer are referred to the relay side (I1, nI2). A single m-section models the 
transmission line. 
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Figure 6. (a) Power transformer single phase equivalent circuit (b) Power transformer current- flux 
characteristics and (c) CTs current- flux characteristics 


3.2. Simulink model of PAD scheme and differential relay 

The developed differential relay with PAD scheme is implemented in Simulink software, 
for evaluating its performance at different operating conditions of the power transformer. A block diagrams 
of the proposed differential relay is shown in Figures 7 and 8. Figure 7 include amper 2pu block to convert 
the phase sampling currents from CT1&CT2 (phase A) to pu-values. The relay +PAD block of phase A 
represents the differential relay (87T) with associated PAD schemes. Scopel in Figure 7 is only used to 
determine the phase angle difference between I1, A and nI2,A. The base current for the HV-side is 5.23 A, 
while the base current for the LV-side is 4.358 A. 
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Figure 7. Simulink Block diagram of protection system (Phase A) 


The model of percentage differential relay elements (phase-A elements) with PAD restraint scheme 
is depicted in Figure 8, where the magnitude values of Idiff (I1-nI2) and 2Ires (11+nI2) are computed using 
‘Fourier analyzer’ Simulink model. The differential relay characteristics given by (3) are represented by 
comparators 2 and 3 along with Simulink embedded function f(u) to generate X and Y signals. 
The differential relay characteristic shown in Figure 1 has boundary values: Iop1 is 0.2 pu , Iop2 is 0.3 pu, 
Ires,min is 0.6 pu and the slop k equal to 0.2. The differential element also involves instantaneous trip 
function using comparator 1 that creates the unrestraint signal UR if the condition: Idiff > 20 pu, is satisfied. 
The PAD scheme is characterized by comparator 4 to produce blocking or releasing signal PAD according to 
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the tripping criterion of the scheme. The function fl(u)= Idiff -0.0001 is used only as a correction function to 
correct some mismatch Simulink values. The output signal of AND1 gate is used to generate 87T differential 
relay signal. Finally, the output of OR gate is used to generate the final trip signal (Trip_A) 
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Figure 8. Differential relay (87T) with PAD scheme 


4. SIMULATION RESULTS AND ANALYSIS 

The power transformer in Simulink model is simulated in normal loading, internal fault and external 
fault as well as magnetizing inrush operating modes. The following Figures 9-13 display the phase A 
secondary currents of CT1 (a), and CT2 (b) in pu values. The same figures show the output signal of 87T 
differential relay (c), PAD scheme (d), and final trip signal of phase A (e). The simulation time of our model 
is 200 ms. 


4.1. Power transformer with normal loading 

In normal loading mode, the connected switches with power transformer are in ‘off? position, 
and the breakers BA and BB are switched on 25 ms after the running program. In this mode of operation, 
the input currents to primary have the same direction of the current flow from the secondary of power 
transformer. For this reason, the output signals of 87T and PAD are zero causing the final trip to be zero, 
as shown in Figure 9. 
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Figure 9. Response of the PAD-based algorithm to normal loading 


4.2. Power transformer with magnetizing inrush current 

The transformer model is provided with a non-linear inductor Lm connected as shown in the 
Figure 6(a) for each phase. Lm is used to account magnetizing current representation. The power transformer 
characteristics (current — flux) of Lm are depicted in Figure 6(b). The magnetizing inrush current can be 
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simulated by switching off all switches and breakers on the LV side of power transformer, and closing circuit 
breaker BA. Figures 10(a) and 10(b) compares the proposed Simulink model and EMTP-RV simulation 
package [1] in the case of inrush currents, where the results were found close to each other. In the obtained 
results of Figure 10(a), the breaker BA was switched on after one cycle of running model. A magnetizing 
inrush current of amplitude up to 8 times the transformer rated current is simulated. It is obvious that such 
inrush current causes the mal-operation of the differential relay. Nevertheless, the output signal of PAD 
scheme is zero logic and the PAD decision is stable and independent of the magnitude of the inrush current. 
Finally, the transformer breakers will not trip. 
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Figure 10. Response of the PAD-based algorithm to magnetizing inrush current, (a) Proposed Simulink 
model, (b) EMTP-RV simulation package [1] 


4.3. Power transformer with turn to ground fault 

Figure 11 presents the performance of the percentage differential relay 87T and the PAD scheme to 
a turn-to-ground fault occurred when the switches Sı &So; are closed. Figure 11(a) shows the results of the 
proposed Simulink model whereas Figure 11(b) illustrates the results of using the EMTP-RV simulation 
package and applying PAD [1]. The results of both methods were similar to each other. The fault is initiated 
at time t = 25 ms with zero fault resistance. Both differential relay and PAD scheme have the ability to detect 
the fault, but with different time instants as appeared in the figure. The different time instants is due to the 
utilization of Discrete Fourier Transform (DFT) for calculating the magnitude values of the differential and 
restraint currents. The final trip signal is released only for the faulted phase, leading to the disconnection of 
the power transformer from the total power system. 
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Figure 11. Response of the PAD-based algorithm to turn-to-ground fault. (A) Proposed Simulink model, (B) 
EMTP-RYV simulation package [1] 
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4.4. Power transformer with turn to turn fault 

Turn-to-turn fault is tested in Figure 12 when Sı of Phase-A is closed at time t = 25 ms. In the 
simulation result both differential relay and PAD scheme can see the fault, but with different time instants. 
In this case, PAD scheme signal is the same final tripping signal and it can lead to disconnect the power 
transformer from total power system with fast response. 
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Figure 12. Response of the PAD-based algorithm to turn-to-turn fault 


4.5. External fault (ACG fault) 

The PAD based algorithm is also tested for external fault (ACG). This fault is located on the low- 
voltage terminals of the power transformer, and it is initiated at 25 ms, with fault resistance 0.001 Q. 
The status of the proposed relay is shown in Figure 13. As seen, the final trip signal is zero logic because the 
CT, and CT? currents flow in the same direction and less than 20 pu. 
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Figure 13. Response of the PAD-based algorithm to AC-G external fault 


5. CONCLUSION 

In this paper the transformer differential protection with PAD based algorithm is modeled by 
Simulink package. In Simulink model, the power transformer and all current transformers are represented by 
real equivalent circuit. Also, non-linear inductance is used to represent the magnetization effects. This 
nonlinearity causes mal-operation in percentage differential relay. The algorithm of transformer differential 
protection relay was developed by adding PAD based algorithm; this improvement is used to solve the mal- 
operation of transformer differential protection on magnetizing inrush currents. The Simulink model is 
simulated at different modes: normal loading, inrush currents and internal and external faults. The simulation 
results show that the proposed algorithm has good reliability, largely independent of harmonic contents in the 


Int J Pow Elec & Dri Syst, Vol. 11, No. 2, June 2020: 1088 — 1098 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 O 1097 


differential current, and transformer parameters. Simulink model shows that the proposed algorithm does not 
require complex computation and can be easily incorporated into existing digital differential relays. 





APPENDIX 
Table 1 Power System Parameters 
Parameters value 
Sourcel: Es =138kV L 0o /50Hz ,R+=7.1Q, L+=53.99mH 
RO=7.1596Q, LO=115.45mH 
Source2: Er =13.8kVL 200 /50Hz ,R+=1.4Q, L+=5.6mH 
RO=1.498Q, LO=11.975mH 
Power 25MVA, 138kV/13.8kV,R1=0.908Q, R2=0.0091Q 
Transformer L1=78.51mH, L2=0.7851mH, Rc=1.19MQ 
Lm=flux/current (pu) / Figure 6(b). 
25VA, 0.25/5 V , 50Hz, R1=2.5uQ, R2=1mQ 
CT1: 100/5 A L1=0.31831uH, L2=0.12732mH, Rc=0.25Q 
Lm=flux/current (pu) / Figure 6(c). 
25VA, 0.020833/5 V , 50Hz, R1= 17.361nQ, R2=1 mQ, 
CT2: 1200/5A L1=2.2105nH, L2=0.12732mH 
Rce=1.7361mQ, Lm=flux/current (pu) / 
Figure 6 (c). 
Line parameters 
RO = 0.1437Q, R1 = 0.3101 Q LO = 11.45 mH, 
L1 = 2.41 mH, CO = 5.635 nF, C1 = 26.8 nF , length= 1km. 
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